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C — ■ Abstract 

■ To study "physical" gauges such as the Coulomb, light-cone, axial or tempo- 
ral gauge, we consider "interpolating" gauges which interpolate linearly between 

' a covariant gauge, such as the Feynman or Landau gauge, and a physical gauge. 

pi 1 Lorentz breaking by the gauge-fixing term of interpolating gauges is controlled 

by extending the BRST method to include not only the local gauge group, but 
also the global Lorentz group. We enumerate the possible divergences of inter- 
polating gauges, and show that they are renormalizable, and we show that the 
expectation value of physical observables is the same as in a covariant gauge. In 
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the second part of the article we study the Coulomb-gauge as the singular limit of 
the Landau-Coulomb interpolating gauge. We find that unrenormalized and renor- 
malized correlation functions are finite in this limit. We also find that there are 
finite two-loop diagrams of "unphysical" particles that are not present in formal 
canonical quantization in the Coulomb gauge. We verify that in the same limit, the 
Gauss-BRST Ward identity holds, which is the functional analog of the operator 
statement that a BRST transformation is generated by the Gauss-BRST charge. 
As a consequence, gA$ is invariant under renormalization, whereas in a covariant 
gauge, no component of the gluon field has this property. 
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1 Introduction 



Although different gauges are formally equivalent, some are simpler than others, or may 
have attractive properties. Covariant gauges are well adapted to perturbative expansion 
and renormalization. However in QCD we are interested in confinement and eventually 
in bound-state problems which are inherently non-perturbative. For such puposes, non- 
covariant gauges such as the Coulomb gauge, the Weyl gauge, the axial gauge or the 
light-front gauge may be attractive. These gauges are considered "physical" in the sense 
that the space of states is believed to be unitary and does not involve ghosts. (For a 
discussion of various gauges see 0.) Indeed non-covariant gauges such as the Coulomb 
gauge and the light-front gauge have recently been used to investigate confinement in 
QCD |, ?]. 

However it is a fact that at the level of quantum field theory, the well established, 
renormalizable gauges for QCD are, on the one hand, covariant, and on the other, in- 
volve "unphysical" particles. These are the fermi-ghosts that are needed to cancel the 
unphysical gluon degrees of freedom. One would like to know whether or not the physical 
gauges really exist in the sense of perturbatively renormalizable quantum field theories, 
and whether they are really unitary in the sense that they may be expressed without 
ghosts, in terms of the two transverse degrees of freedom of the gluon. We shall see that 
for the Coulomb gauge, the answer to the first question is "yes" and to the second, a 
slightly qualified "no". 

The point of view which we adopt in the present article is that the BRST formulation 
provides a reliable method of quantizing and perturbatively renormalizing non-Abelian 
gauge theories. (For a review see |4j and ||.) The existence and properties of physical or 
canonical gauges will be investigated deductively starting from the BRST formulation. 
To be sure, this inverts the historical order in which gauge theories were first canonically 
quantized, and subsequently the BRST method was found; however the canonical method 
has remained heuristic, and to this day does not allow systematic renormalization. 

There are two different problems raised by the commonly used "physical gauges": 
(i) the breaking of Lorentz covariance and (ii) an arbitrariness due to incomplete gauge 
fixing. For example the Coulomb gauge condition V • A = obviously breaks Lorentz in- 
variance. It is also an incomplete gauge-fixing in the sense that it leaves a one-parameter 
family of gauge transformations arbitrary, namely gauge transformations g(t) that are 
independent of the spatial coordinate x, but may depend on the time t. Similarly, the 
Weyl gauge condition A = leaves arbitrary a 3-parameter family of gauge transforma- 
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tions, g{x). We call the dimension a of this parameter space the "degree of arbitrariness" 
of the gauge, and we have a = 1 for the Coulomb gauge, and a = 3 for the Weyl, the ax- 
ial, and the light-front gauges. Not surprisingly, the degree of arbitrariness of the gauge 
determines the dimension of the divergences of Feynman integrals that are not controlled 
by usual ultraviolet regularization. 

Strictly speaking, incomplete gauge fixing with a > implies that the correlation 
functions of charged fields actually vanish at generic space-time separation. For example 
in the Coulomb gauge, the arbitrariness under g(t) implies that the correlation function 
of two charged fields vanishes at unequal times, {i/j(x, t)ip*(0, 0)) = for t ^ 0, even in 
abelian gauge theory. This vanishing of correlation functions due to gauge arbitrariness is 
not what one has in mind by a 'physical' gauge, and it is usually overcome in continuum 
gauge theory by additional gauge fixing by more or less explicit prescriptions.^ Incom- 
plete gauge fixing would appear to be the origin of ambiguities that occur in higher loop 
diagrams 0, and which make the formal Coulomb gauge, defined by canonical quanti- 
zation after elimination of the Coulomb-gauge constraints, not particularly well-defined. 
Consequently it is very misleading to speak of the Coulomb gauge, as in the question, 
'What are the Feynman rules for the Coulomb gauge?'. Unless one is willing to accept 
the vanishing of correlation functions of charged fields at unequal times, this question 
cannot have a unique answer without further stipulation of the gauge condition. This 
applies to all gauges with a > 0. 

We shall deal with both problems, Lorentz breaking and gauge arbitrariness, by the 
device of an "interpolating" gauge. For example the gauge condition — adoA + V- A = 0, 
with < a < 1 interpolates between the Landau gauge, a = 1, and the Coulomb 
gauge, a = 0. For a > the gauge condition is regular, in the sense that the degree of 
arbitrariness vanishes, a = 0, but Lorentz invariance is broken for a ^ 1. This allows 
one to first address the problem of Lorentz breaking in a regular gauge, and then to see 
if the singular limit a — > yields finite correlation functions. In the present article we 
shall use and extended BRST symmetry to control the violation of Lorentz invariance, 
and we shall then study the Coulomb gauge limit of the Landau-Coulomb interpolating 
gauge. 

Use of an interpolating gauge and an extension of BRST symmetry to control the 
violation of Lorentz invariance, was reviewed by Piguet J7|, particularly for the interpo- 
lating light-cone gauge. Doust || used a gauge which interpolates between the Coulomb 
and Feynman gauge to regularize the Coulomb gauge, and showed that extra terms in 

■"■In lattice gauge theory, gauge-fixing is frequently left incomplete. 
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the Feynman rules which he obtained in the Coulomb-gauge limit correspond to an ad- 
ditional potential term obtained by Christ and Lee from an operator ordering of their 
Coulomb Hamiltonian [ID]. Difficulties with renormalization in the Coulomb gauge were 



exhibited by Doust and Taylor 0. The Weyl gauge (A = 0) has been studied by Rossi 
and Testa ||11|| , and by Cheng and Tsai [0 . 



As commonly used in non-Abelian gauge theories, BRST-invariance provides a substi- 
tute for invariance under local gauge transformations which is broken by the gauge-fixing 
term. In Lorentz-covariant gauges, one uses the BRST method to enumerate the inde- 
pendent divergent counter-terms necessary to ensure finitness of the renormalized theory, 
while preserving all requirments of gauge invariance for physical quantities. It is a power- 
ful algebraic method of great generality, relying as it does on the simplicity of invariance 
under a generator s that is nil-potent s 2 = 0. 

In the first part of this article, we develop an extension of the BRST method that also 
provides a substitute for invariance under global Lorentz rotations when the gauge-fixing 
term breaks global Lorentz invariance as well as local gauge invariance. The method is of 
considerable generality in that it does not rely on particular properties of the symmetry 
which is broken by the gauge-fixing term, but only that the symmetry operations form a 
Lie group, and it allows us to explicitly enumerate all counter-terms. 

For the class of interpolating gauges, defined by (ad)^A^ = f, with a a non-singular 
matrix, the partition function is formally given by the Faddeev-Popov formula 

Z = J dA6[(adrA, - f] det[(adrD,(A)] exp[-S YM ], (1) 

where Sym is the Euclidean Yang- Mills action. Feynman graphs contain denominators 
of the form k^a^ky and k 2 . As long as a is non-singular, these denominators provide 
the same degree of convergence in all directions in fc-space as the corresponding denom- 
inator k 2 in covariant gauges. Consequently in this class of interpolating gauges, power 
counting of graphs is exactly the same as in Lorentz-covariant gauges. The problem 
of renormalizability is reduced to an algebraic one of enumerating the form of possible 
local divergent terms, which we control by extended BRST-invariance. On the contrary, 
because of gauge arbitrariness in the limiting cases of the Coulomb, light-cone or other 
singular gauges, the degree of convergence depends on the direction in k-space, and a 
more detailed analysis is required to determine if the limit is finite. 

In the second part of the present article, we analyse the singular Coulomb gauge limit, 
a —>■ 0, from the Landau-Coulomb interpolating gauge. For this purpose we express the 
partition function Z as a functional integral in phase space, and then make a linear shift 
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in the field variables in order to exhibit a symmetry (r-symmetry) between the fermi 
and bose unphysical degrees of freedom. Individual closed fermi-ghost loops and closed 
unphysical bose loops diverge like a -1 / 2 , but they cancel pairwise by virtue of the r- 
symmetry. Consequently the correlation functions are finite in the limit a — > from the 
Landau- Coulomb gauge. This remains true for the renormalized correlation functions. 
(See remark 1 at the end of sect. 9.) However we also find that there are one- loop graphs 
that vanish like a 1//2 , and that are missing in the formal (a = 0) Coulomb gauge, but 
which cannot be neglected because they give a finite contribution when inserted into the 
graphs that diverge like a -1 / 2 . It remains a logical possibility that these two-loop graphs, 
that are missing in the formal Coulomb gauge, are mere gauge artifacts that decouple 
from expectation values of all gauge-invariant quantities such as a Wilson loop. However 
there is at the moment no argument to show that they do. Indeed unless for some reason 
these two-loop graphs decouple from all physical amplitudes, then the ghosts do not 
decouple in the Coulomb gauge limit, and the Coulomb gauge is not unitary in the usual 
sense of being a canonical theory of the transverse gluon degrees of freedom. 

Nevertheless we find that correlation functions of the Coulomb-gauge limit of the 
Landau- Coulomb interpolating gauge do exist, and moreover they display a kind of sim- 
plicity that is absent from covariant gauges. A certain Gauss-BRST Ward identity holds 
in the Coulomb gauge limit which implies, among other things, that the time-time compo- 
nant of the gluon propagator g 2 Doo is a renormalization-group invariant and thus depends 
only on a physical mass, Aqcd, but not on the ultra-violet cut-off, A, nor the renormal- 
ization mass, fi, which may make it a useful order parameter for color confinement. No 
component of the propagator has this property in a covariant gauge. 

2 Interpolating Gauges 

In this section we introduce interpolating gauges for various familiar classical gauges. 
The Landau and Coulomb gauges are defined by — d A + V ■ A = and V ■ A = 0. The 
Weyl and axial gauges are frequently defined by A = and A 3 = respectively. However 
if periodic boundary conditions are introduced in time or space, the conditions A = 
and A 3 = are too strong, and cannot be maintained. For they fix to unity the values 
of straight-line Wilson loops trPexp(J dx^A^) that close by periodicity, which however 
are gauge-invariant objects. We take instead as the Weyl and axial gauge conditions the 
weaker conditions d A = and d 3 A 3 = 0. In momentum space these read k A (k) = 
and k 3 A 3 {k) = 0, so the weaker conditions differ from the stronger ones by zero modes 
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only. Similarly for the light-front gauge condition, instead of — Aq + A3 = we take 
(d + d 3 )(-A + A 3 ) = 0. 

All these gauge conditions have the linear form (Pd) ■ A = 0, where in the various 
cases P^ is the projector 

Landau : P^ = 5^ = diag(l, 1, 1, 1) 
Coulomb : P; = diag(0, 1, 1, 1) 
Weyl:P;=diag(l, 0,0,0) 
axial : P^ = diag(0, 0, 0, 1) 
light - front : Pf = P 3 ° = P 3 = P ° = 1/2 and P; = otherwise. (2) 

These projectors have a null space of dimension a = 0, 1, 3, 3, and 3 respectively, where 
a is the degree of arbitrariness of the gauge, as defined in the Introduction. 

To separate the problem of violation of Lorentz invariance by the gauge-fixing con- 
dition from the problem of the arbitrariness of the classical gauges, we introduce an 
interpolating gauge defined by the condition (ad) ■ A = 0. Here a is the numerical 
matrix 

a = P + aQ, (3) 

where P is one of the above projectors, Q = (1 — P) is the orthogonal projector, and a 
is real, in the interval < a < 1. These gauges interpolate between the Landau gauge, 
at a = 1, and any one of the above singular classical gauges, which is achieved at a = 0. 

For the quantum field theory we consider the slightly more general gauge condition 
(ad) ■ A = f. By the usual Faddeev-Popov argument, the partition function, eq. ([1]), is 
expressed in terms of the local Faddeev-Popov action, 

S FP (A, c, c) = S YM (A) + Jd*x{ (2(3)~ 1 [(ad) ■ A] 2 + (ad)c ■ D(A)c } (4) 

where D(A) is the gauge- covariant derivative [D^(A)c\ a = d^c a + f^A^ c d , and (3 is a 
gauge parameter. 

From this action, one reads off the ghost propagator 

G = -i (k-k')~\ (5) 

where 

k' = ak = Pk + aQk. (6) 
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Similarly the gluon propagator D is obtained from the quadratic part of the gluon action 
(1/2)(A,KA) by K^D^ = -i8*. From the Faddeev-Popov action we have 

K iw = k 2gv> _ k » k v + /r 1 W, ( 7 ) 
and one easily verifies that the gluon propagator is given by 

= -% k- 2 [ 9fiu - (k ■ k'Y\kX + k'^ky) + (k ■ k')~ 2 ((3k 2 + k' 2 )kM . (8) 

As long as a is a non-singular matrix, namely for a > 0, convergence of Feynman integrals 
is independent of direction in momentum space. The familiar power counting arguments 
hold, and Feynman integrals may be regularized by dimensional regularization. 

We now consider some special cases. A Landau-type interpolating gauge is obtained 
at /3 — 0. In this case the propagator satisfies the generalized transversality condition 
k'^D^y = 0. For a = 1, we have the Landau-gauge propagator, so this gauge interpolates 
smoothly between the Landau gauge and the classical singular gauges. 

A Feynman-'tHooft type gauge is obtained by choosing (3 so that the double pole 
becomes a simple pole. For Coulomb, Weyl and axial gauges, the projector P commutes 
with the metric tensor g = diag( -1, 1, 1, 1), and we have (Pk) ■ (Qk) = 0. In these 
gauges the double pole is eliminated by setting j3 — a, for we have 

(3k 2 + k' 2 = a[(Pk) 2 + (Qk) 2 } + (Pk) 2 + a 2 (Qk) 2 = (1 + a)k ■ k', (9) 

which gives the propagator 

= -i k- 2 { 9llu + (k ■ k'Y\-kX - k'^K + (1 + a)k^k u }} . (10) 

This propagator has the attractive 'tHooft-type property that it is block diagonal in the 
P-Q subspaces. 

On the other hand, for the light-front gauge (Pk) ■ (Qk) ^ 0, but Pk is a null vector, 
(Pk) 2 = 0. In this case the double pole is eliminated by setting (3 = a 2 , for we have 

pk 2 + k' 2 = a 2 [2(Pk) ■ (Qk) + (Qk) 2 } + 2a(Pk) ■ (Qk) + a 2 (Qk) 2 

= 2ak-k', (11) 

which gives the propagator 

= -% k- 2 { 9llu + (k ■ k'Y\-kX - k'k u + 2ak lx k v }} . (12) 
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In the last two expressions for D^, the Feynman gauge is obtained at a = 1, so 
these gauges interpolate smoothly between the Feynman gauge and the classical singular 
gauges. 

We write these expressions explicitly for interpolating Coulomb gauges. In this case 
we have k' = Pk + aQk = (ak , k), and the ghost propagator is given by 

G=-i ^— — - . (13) 

For the gauge which interpolates between the Landau and the Coulomb gauges, the gluon 
propagator is given by 

_ J_ / kjkj\ kjkj a 2 kl 
13 ~ k 2 [ 9 » P ) P 02 _ ak 2 ? ^ 

= _ ak k t 

{k 2 -aklf V ' 

p 

i D 00 = — . (16) 

(fc 2 - afcg) 2 

These expressions are obtained by partial fractionation, and there is no singularity at 
k = for a > 0. 

It is easy to understand intuitively how the Coulomb-gauge limit from the Landau- 
Coulomb gauge fixes the gauge arbitrariness of the Coulomb gauge discussed in the Intro- 
duction. Under the residual gauge freedom of the Coulomb gauge, A transforms accord- 
ing to Aq — > g^(t)A g(t) + g^(t)dog(t), where the inhomogeneous term is x- independent. 
With periodic boundary conditions, the Landau-Coulomb gauge condition adoA = V- A 
for a > gives d J d 3 xA = 0. However, as one sees from the above expression for the 

— * 

A — A propagator, D 00 vanishes at k = for all finite a, so the stronger condition 
/ d 3 xA = in fact holds in the Landau-Coulomb gauge for all finite a. This provides 
the additional gauge-fixing condition needed to make the limit a — > well defined. By 

— * 

contrast D Q0 in the Feynman- Coulomb gauge, given below, becomes ill defined at k — 
with periodic boundary conditions, in the limit a — > 0. 

For the gauge which interpolates between the Feynman and the Coulomb gauges one 

has 

(17) 




akl 



(18) 
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kp 1 Qik^ 

There is no mixing of space and time components of the gluon propagator in this gauge. 

These expressions for the propagators are quite illuminating. The transverse part of 
Dij is the Coulomb gauge propagator. The parameter a acts as a regulator for simul- 
taneity in the Coulomb gauge. These expressions imply exact compensations between 
the "unphysical" contributions in internal loops between gluon and ghost propagators. 
The main thing is of course that these compensations occur because the poles of the 
propagators of the unphysical fields with opposite statistics sit at the same point, i.e, at 
P - akl = 0. 

For completeness, we indicate the form of propagators in the interpolating gauges for 
the light-front gauge quoted in ||, with k'^ = (a^pck)^ given by 

k^ = ak u k' 2 = ak 2 , k' 3 = ^-^k 3 + ^-^k , k' = ^-^k 3 + ^-^k . (20) 

We have k ■ k' = a{k\ + k 2 ) + (fcf — kfy. For the gauge which interpolates between 
the light-front and Feynman gauges, namely, with /3 = a 2 which eliminates the double 
pole in the gluon propagator, the above expression for the gluon propagator reads 

5 a 

D = -i — 

3 k 2 



A- = 

D i+ = i (1 - 
J D__ = 

D - — 
D+ -~ k-k> 



kj{h + fc ) 
k 2 k-k' 



. (ki + k Q ) 2 , , 

D ++ = 2i(l-a) { ^ k ^ , (21) 

where i,j = 1, 2, D^± = D^ 3 ± D^q, and D±± = D 3 ± ± D ±. We observe that _D M _ = 
at a = 0, and the light-front gauge condition is satisfied. 



3 BRST symmetry for local gauge and global Lorentz 
invariance 

Suppose that we have a Lie algebra with basis Xi and structure constants fijk, so 
[Xi, Xj] = fijkXk. According to the BRST method, for each generator Xi we intro- 
duce a corresponding Grassmann or ghost variable q. The BRST operator s acts on 
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these variables according to 

sCi = —-fijkCjCk- (22) 

It is nilpotent, s 2 = 0. The preceding relation is isomorphic to the action of Cartan's 
exterior differential operator d acting on the Maurer-Cartan form u = Uit 1 = g~ 1 dg of 
the Lie group. 

We wish to apply this method to the Lie group which consists of local gauge transfor- 
mations and global Lorentz transformations. The structure constants of this group are 
given by 

[G a (x), G\y)} = f abc 5{x - y)G c {x) 
[H^G a (x)\ = -( Xfl d v - x v d^)G\x) 
[H K \, H^} = gx^H KU — g KfM Hxu — g\uH Kfl + g KU Hx^. (23) 

According to the method described above, corresponding to the local generators G a (x) 
we introduce the usual ant i- commuting Grassmann field variables c a (x) and correspond- 
ing to the and we introduce the global generators V^ u = —V Vfl , so Cj = (c a (x), V^). 
In 4-dimensional space-time there are 6 independent generators V^ v . For the structure 
constants just found, the BRST operator s acts according to 

sc a (x) = -^f abc c b (x)c c (x) + z V^x»d u c a (x) (24) 

stV = -W/, (25) 

where the parameter z will be determined shortly. Because V^„ is a Grassmann variable, 
{V 2 )ku = V R \g Xll Vfj, v is an anti-symmetric matrix, (V 2 )^ u = — (V 2 ) u ^. Equation ^4] deter- 



mines the normalization of the ghost field c a (x), and eq. |25| determines the normalization 
of Grassmann variables V^ v . The parameter z is most easily determined by requiring 
that s be nil-potent, s 2 = 0, which gives z = —1. We could as easily have derived the 
corresponding result for the Poincare group.[| 

The BRST operator associated to the Lie algebra just defined is of the form s = s 9 +sl, 
where s g and sl satisfy (s g ) 2 = (sl) 2 = s 9 sl + sls 9 = 0. On the fields c a (x) and 
they act according to 

Sg c a {x) = -\f abc c b {x)c c {x) 
s L c a (x) = V^x^d u c a {x) 
s g Vx v = 

s L Vx u = ~(V 2 )x U . (26) 

2 Use of a "large" BRST operator in the present context was suggested to us by Massimo Porrati 
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The BRST oprator s defined here may be viewed as a "large" BRST operator, which is 
the usual BRST operator s g for the local gauge group extended by the BRST operator 
sl for the global Lorentz group.f] 

To determine the action of the BRST operator on the connection A )J a (x), we could 
start with the familiar transformation law of the connection under local gauge and global 
Lorentz transformation, 

G(u) A*{x) = d^ a {x) + f<*°A*(x)u%x) 

H(e)A fl a {x) = e K x x K d x A^ a (x) + e^A v a (x). (27) 

A more economical way is to construct the most general expression sA^ a = (s g + s^A^ 0, 
which satisfies s 2 = 0. Suppose that si acts according to 

s L A^ a (x) = Zl V K x x K d x A^(x) + z 2 V fl u A u a (x). (28) 

Here z\ and z 2 are parameters that are determined by the condition (s_l) 2 = 0. From 
eq. (f26|) one obtains z\ + (zi) 2 = and z 2 + (-22) 2 = 0. We take z\ = — 1 because Z\ — 
gives a trivial field transformation law. The components of A^ transform either like 
scalars (z 2 = 0) or a vector [z 2 = —1). We take the vector case and obtain 

s L A, a (x) = -V K x x K d x A^(x) - V/A u a (x). (29) 

Finally, suppose that s g acts on A according to 

s fl V = + Z2/f abc A u b c c . (30) 

We obtain from {s g ) 2 = that z 2 ^ u = 5^ u , and Z\^ remains arbitrary. The condition 
s 9 sl + sis g = gives ziV = Vz\. Because is arbitrary, z\ is of the form z\p = zS^ u . 
We write A^ a = zA' a , and obtain for the BRST operator s = s g + sl, 

sV = <V a + f abc A, b c c - V K x x K d x A„ a (x) - V^A v a {x), (31) 

where we have dropped the prime on A'. This completes the determination of the action 
of the BRST operator on the basic fields A and C. 

3 The Lorentz rotations are a subset of general reparametrization. If we define the vector £^ = VJfx 1 ', 
out of the constant ghosts V and the coordinates x^, we have = £, p d p ^ and slc = ^d^c. 
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4 Extended action 



The partition function in eq. (]T|) may be expressed in terms of the local Faddeev-Popov 
action, 

S FP ($) = S YM (A) + J d 4 x[-(adybA, + (adYcD(A) tl c+^b 2 }, (32) 

where D(A) is the gauge-covariant derivative [_D At (A)c] a ee d^c a + f^A^ c°, and $ rep- 
resents the set of fields $ = (A, c, c, b). We introduce a corresponding set of sources, 
J — (J a, J ci J ci Jb)i an d write 

Z(J) = J d$exp[-S F p($) + ($, J)], (33) 

where d$ = dAdcdcdb, and 

($, J) ee J d A x(A ■ J A + c- J c + c- J 5 + b- J b ). (34) 

The Faddeev-Popov action is not invariant under Lorentz transformations because of 
the appearance of the numerical matrix a^ v '. Consider instead the extended action 

S ext ($, V) ee S YM (A) - s I d A x[(adrcA^ - ^cb], (35) 

where s is the "large" BRST operator that expresses the substitute gauge and Lorentz 
transformations. Its action on A, c and V is defined in eqs. (|3T|), ( |24] ) and (|25|), and its 
action on c and 6 is defined by sc = b and sb = 0, which preserves s 2 = 0. Because 
the Yang-Mills action Sym{A) is both gauge and Lorentz invariant, it is invariant under 
the "large" BRST operator sSym(A) = 0, and consequently so is the extended action, 
sS ex t{&, V) = 0. The extended action differs from the Faddeev-Popov action by terms 
linear in the global Grassmann variables V introduced in the preceding section, 

S ext ($, V) = S FP ($) - J d^xiady-ciV^x^dxA, + VJA V ). (36) 

We treat the variables as external sources,^ and define the extended partition function 

Z(J,V) ee J rf$exp[-5 cxt ($,K) + ($, J)}. (37) 

The original partition function is obtained from it by Z(J) = Z(J, 0). Because there are 
6 independent global Grassmann variables V, there are, in all, 2 6 terms in the expansion 
of Z(J, V) in powers of V. They are related by the symmetry generated by the large s. 

4 Equivalently we may treat the vector = V^x" as an extended source. 
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The usual argument that the expectation values of gauge-invariant observables are 
independent of the gauge parameters must be slightly modified because the variable V 
is not integrated over. We consider only s-invariant observables W are indepndent of V. 
We shall show that (W) is independent of the a matrix when the external source V is 
set to 0. We also set all sources J to 0, and we have 

d(W)/doT = jd$W Jd*x s(d u c A M )exp(- 1 S cxt )| y=0 

= Jd A x J rf$ s[W d v c A^exp(-S cxt ) ]\ v=0 , (38) 

where we have used sW = sS cxt = 0. At V = we have s = s g , where s g is a derivative 
with respect to the variables of integration <3> = (A,c,c,b). This gives d(W) / da^ v = 0, 
as asserted. We conclude that for physical observables, the interpolating gauges gives 
the same expectation values as the covariant gauges. In particular they are independent 
of the gauge parameter a, and similarly for j3. 



5 Quantum Effective Action 

To exploit BRST symmetry in renormalization theory, it is helpful to also introduce 
sources for the BRST transforms that are non-linear in the fields. We therefore define 
the (fully extended) action 

£($, V, K, L, M) = S cxt + (K, sA) + (L, sc) + M ■ sV, 

= S cxt + s[-(K, A) + (L,c) + M-V] (39) 

where K fl a (x) and L a (x) are the usual sources for sA fl a (x) and sc a (x), and we have intro- 
duced a corresponding source M* v = —M Ufl for sV = —V 2 , with M ■ sV = \M^ v sV^ u . 
These sources are not acted on by s, sK = sL = sM = 0. The action S is invariant 
under the "large" BRST operator, sS = 0. 

We define the corresponding partition function 

Z(J, V, K, L, M) ee J d$exp [-£ + ($, J)]. (40) 

It satisfies 

m= = - (y2) - z ' (41) 

The BRST operator s has been defined as a linear differential operator that acts 
on (and mixes) the variables $ = (A, c, c, b) and V. Because only the $ variables are 
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integrated over, it is convenient to decompose s according to s = s$ + sy, where s$ acts 
only on the $ variables, and sy only on V, so s$V = sy§ = 0. The explicit form of sy 
is 

By the invariance of £ with respect to s = s$ + sy, we have 

W-§ = IHJa,±)-V~±)-V»±)\z (43) 

The free energy W(J,V, K, L, M) = In Z(J, V, K, L, M), satisfies the corresponding 
equations 

m 

sV^ (44) 



8M^ 

< J * jk ] + ^ tl> + v» + {sV) -w = °- (45) 

We make a Legendre transformation from the variables J = (Ja, J c , Jc, Jb), an d the 
free energy W to the external field variables $ = (A,c,c,b), and the quantum effective 
action T, 

r($, V, K, L, M) = ($, J) - W(J, V, K, L, M) (46) 

where 

5W 5W 5W , SW 

= c = c = o = 

5Ja„ SJ c 5Jc 5J b 

(47) 



sr sr sr 

1 ~ SA„ Jc ~ 5c c ~ 5c b ~ 5b 



(48) 



and 



5T 
5V 



5W 



sv 



5T 
5K 



5W 



5T 
5L 



5W ST 



5W 



5K 1 5L 5L 1 5M 5M 1 ^ 
Here and elsewhere, all derivatives with respect to fermionic variables are left derivatives. 
In terms of T, eqs. (44) and fl4"5|) give 



5T 

SM^ 



sV, 



<V 2 



(50) 
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and 

.sr 5T. .sr sr. ,5r sr 5r n 



This type of equation, which was introduced in |13|, now includes a V-M term. Here 



is it assumed that there is no gauge anomaly. No Lorentz anomaly can occur in D = 4 
dimensions. 

Because the gauge condition is linear, we may solve the equations of motion to obtain 
the dependence of T on the Lagrange multiplier fields c and b. As this is standard, we 
simply give the result WQ, 



r(A, K, c, L, c, b, V, M) — J d A x[-(adYbA^ + {^)b 2 } + f(A, K + (ad)c, c, L, V, M) (52) 

The property that the K and c dependences are only through the combination K^ + ad^c 
can be imposed as a Ward identity in the class of linear gauges that we consider. This 
plays an important role in the renormalisation program. 

The master equation satisfied by r(A, K, c, L, V, M) is symmetric in the pair V, M 
and the other variables, 

( M' JK ] + ( ' Jl ] + 6V ' 6M ~ °- (53) 
T has the simple dependence on M given by 

df 



dM^ v 



sV^ = -(VV- (54) 



6 Form of Divergences 

The new V — M term has the same structure as the other terms, so we may use familiar 
arguments, which we now sketch, to determine the form of possible divergences to each 
order in h, when using a regulator that preserves Lorentz and gauge symmetries. We 
make a loop or h expansion of T, using any suitable regularization for divergences. 

f = ^F\ (55) 

n 

To find r°, we observe that E, eq. (|39|) , is of the form 

E(A, K, c, L, c, b, V) = J d 4 x[-(ad)"bA^ + ^f3b 2 ] + E(A, K + (ad)c, c, L, V, M) (56) 
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where 

E(A, K, c, L, V) = S YM + (K, sA) + (L, sc) + M ■ sV. (57) 

This gives 

f° = E. (58) 

We will impose that V is renormalized while satisfying the master equation (|53]). From 
the s-invariance of SVm ; and from sA = 5Y>/5K, sc = 5Y>/5L, and sV = 5H/5M, we 
have 

( M' Jk ) + [ Jc~'Jl ) + WJm~" (59) 

so the master equation (^3|) is satisfied by T = E. We define the star product 
f „f - fST* Sf b 5T a 5T b 5T a ST b 

To each order n in ft, eq. (E3I) reads 



^ f p *f 9 = 0. (61) 

p+q=n 

We assume that renormalization has been done to order n — 1, so that T p for p = 
0, • • • (n — 1) is finite, and that eq. (|61~D is satisfied to order n — 1. We separate the regular 
and divergent parts of order n, 

f » = f I + f » v , (62) 

where the first term is the renormalized part of the n-th order effective action and is 
finite. By hypothesis, the only divergence in eq. fl6~T|) comes from T^ iv - The divergent 
part must satisfy eq. (|6l|) separately, namely 

*f div = 0, (63) 
where the linear operator a, defined by 

aT = t*T + T*t, (64) 

has the explicit expression 



, j4 /(5S 5 5E 5 <5E 5 5E 5 8E 5 5 

a = (TX\ —— + 77777 + 7T7-7- + "7" 7^ + 7TT7777 + 



5c 6c 5L 5M6V 5V 5M 



(65) 
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It is nilpotent a 2 = 0. Here a represents the symmetry of the "large" BRST operator, 
with the obvious decomposition into local gauge and global Lorentz parts, a = a g + cr^, 
that corresponds to s = s g + sl- From eq. (Q) we have 

— = 0; „>1, (66) 

so T ra is independent of M. 

Consistent with the last equation, with locality of divergent terms, with global color 
invariance, with the ghost quantum numbers (0, —1, 1, —2, 1) and dimensions (1, 2, 1, 2, 1) 



of the variables (A^, K a '^^ c a ,L a } V^ u ) on which T n depends, eq. (|63|) has the solution 

(67) 



rSiv = / d*x 



Cl -Xu + ^(K^c 2 /Al + c,L a C a ) 



4 



where ci, c<i^ and C3 are divergent constants of order H n . The operator x^dy — x^d^ 
may appear in T n in the combination V^x^dy. However x^d v — x v d^ is dimensionless 
and carries no ghost or global color quantum number so it does not affect our counting 
arguments, which exclude the explicit appearance of V in the last equation. However a 
V dependence is introduced into r n from the definition of a, so that V appears in the 
expansion of the a-exact term. 

With this result we have achieved our goal of limiting the number of possible diver- 
gences, by maintaining invariance under the larger group of substitute gauge and Lorentz 
invariance. Indeed only the combination / d 4 xF 2 u is invariant under a without being ex- 
act, of the form aX. (J d xF 2 is said to be the cohomology of the operator a.) For if 
only invariance under s g or a g were enforced, then the most general cohomology would 
be / d 4 x(cEE 2 + cbB 2 ) , if ordinary rotational invariance is preserved by the gauge fixing, 
where ce and cb are independent renormalization constants. Indeed, E 2 and B 2 are sep- 
arately invariant under s g , and in 0, it was necessary to assume ce = cb- This is now 
established for the gauges considered here. On the other hand the breaking of Lorentz 
invariance by the gauge fixing does lead to the Lorentz non-invariant divergent terms 
a K a '^ 'A a „, which however are exact cr-forms. 

If ordinary rotational invariance is maintained by the gauge fixing, then C2, fl u is a 
diagonal tensor with C2,\ l = c 2> 2 2 = c 2> 3 3 7^ c 2> 4 4 . In Lorentz-type gauges, defined by 
(3 = in eqs. ([32]) or fl35D, the (possibly) divergent constant c 3 vanishes, c 3 = 0, by 
virtue of the factorization of the external ghost momentum, as it does in the Landau 
gauge ITS . 
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7 Multiplicative Renormalization 



In previous sections we had implicitly absorbed the coupling constant g into the SU(N) 
structure constant f abc . Since we are interested in the perturbative expansion we now 
make the coupling constant explicit by the substitution f abc — > gf abc . 

We define S r ($, V) r as the local part of V = V R + T div . We call E r ($, K$, V) the 
renormalized action since by inserting expE r ($, K$, V) in the path integral over the $, 
all relevant local divergent counterterms are present to determine finite Green functions 
of the fields A, c, c, b and of their BRST transformations which satisfy the BRST master 
equations. 

The result found in eq. fl6?D, proves that the renormalized action S r ($, V) r has 
the following form: 

E r ($,lf*,V% = J d 4 x{ i| F, u (Z^A u ,Z g g) | 2 

+ Z c [ {adYc + K"]Ztf( d v c + Z g g[ Z P A ^A P , c] ) 

- \z c Z g g\ c ,c]L + ^b 2 + b{adYA» - \m^(V 2 ), v } (68) 

For the sake of notational simplicity, we use the graded commutator notation, [X, Y] a = 
f b a c X b Y c , and F^A, g) = d IM A u - d v A„ + g[A^ A v \. 

The relation between the renormalization constants Z and the constants c appearing 
in eq. (pTj) is 

Z c = l + c 3 + | 

^ = 1 -| ( 69 ) 

Thus, the effect of renormalisation, constrained by the BRST invariance, can be seen as 
the following redefinitions of fields and parameters: 

A^ — > Z v Ail A v 

c— > Z c c 

9^Z g g 
K" -> Z A ^K V 

Z c 
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b^b 

M-^M (70) 

Here and transform contragrediently under renormalization, as do c and L, so 
that the master equation is invariant under renormalization in any finite order. 

Equation ( [FTP shows that the renormalization is (matricially) multiplicative for the 
fields, sources and parameters of the theory and that, as compared to the covariant case, 
the breaking of Lorentz invariance by the gauge fixing term induces a mixing by the 
renormalization of the 4 components in and K^. Let us stress that the simplicity 
of the renormalization of c b, (5 and which generalizes that of covariant renormaliz- 
able gauges, is a particularity of linear gauges for which one can maintain the K and c 
dependences through the combination K + adc. 

These equations indicate the existence of a renormalized BRST symmetry for the 
action E r ($, K$,V), in eq. (f||). We will shortly display its expression. It is however 
instructive to rederive the renormalized action S r ($, K§, V), using the method displayed 
in which has the advantage of determining at the same time the renormalized BRST 
invariance of the theory. 

In this method, one parametrizes the renormalized action E r , including all relevant 
couterterms, as 

E r = £ r ($,V) s r $)+Ms r V (71) 

Recall that <£> stand for all fields, A, c, c, b. One has assumed that the dependence on 
the sources K's of the BRST transformations is linear, which will be checked by self 
consistency. Then, s r $ stand for field polynomials in the fields $, which can be expressed 
as the action on $ of a yet undetermined graded differential operator s r . 

One can show that the content of the Ward identities of the BRST symmetry is that 
(i) S T is invariant under the action of s r and (ii) s T is a nilpotent operator 

s 2 r = (72) 

and 

s r S r = (73) 

To compute the possible action of s r $ , with 5^$ = 0, one uses the results of sec- 
tion (3). Up to inessential overall factors, the only freedom left in determining the action 
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of s r is a matricial redefinition of A^, that is, A^ — > Z^A U , and the rescaling g — > Z 9 g so 
that the requirement s 2 = implies 

SrC=-^f[c,c]-V»X»d u C 

Si z;a v = d,c + z g g\z;A u , c] - v£x x d K z;A u - v?z;a v 

(74) 

and 

s r c = Z~ x b 

s r b = 0. (75) 

This allows one to identify Z\ v = Z v . Notice also the freedom in rescaling the field b. In 
the expression of the action, one furthermore sees that a rescaling of b only amounts to 
a rescaling of the partition function Z, which is unobservable. 
We also remark parenthetically that one can write 

s T A, = %c + Z g g[A„ c] - V^x' x d' K A, - V; A, (76) 

where d'^ = Z~ lv d U) x' M = Z^x v and = Z^VJf Z^ lu . This shows another interesting 
property of the class of non-covariant gauges that we have introduced: the transformation 
of the components of A implied by breaking of Lorentz invariance, while maintaining 
BRST invariance, can be absorbed into a transformations of space time coordinates, 
x — > x', together with redefinitions of constant ghost matrix elements V —>■ V . (One has 

Sr y = -v'V.) 

The non-trivial part of the cohomology of s r with dimension 4 is | d^A v v — d v A r ^ + 
Z g [A v ^ A r u ] | 2 , with A^ = Z^A U ; the rest of S T can only be s r -exact terms with dimension 4. 
By using the anti-ghost equation of motion as a Ward identity, which implies that no 
quartic ghost interactions occur in the action, together with the property that the b- 
dependent part of the action does not need counter-terms, one concludes that S T must 
be of the form 

S V = J d'x{ ~| d,Al - d v A^ + Z g [A^, Al] | 2 + 

+Z c s v {c[ (adTA^ + t h ) } ) (77) 

If we now expand S T , using the definition of s r , and insert this into eq. ([HP, we exactly 
recover the formula giving S r ($, K®, V) in eq. (|68|). 
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The renormalized action S r ($, V), which is suitable for the computation of renor- 
malized Green functions, is thus invariant by construction under the action of the oper- 
ation s r , which is the renormalized expression of BRST symmetry. 

8 Gauss-BRST Ward Identity 

In the remainder of the article we shall study the Coulomb gauge limit. The preceding 
results hold in particular for interpolating Coulomb gauges, when the matrix a is diago- 
nal, and «oo = a and an = 1. As compared to the case of covariant gauges, there is just 
one extra renormalization constant, with A4 — > Z^Ai and Aq —>■ Za Aq, and Z ^ 7^ Z Ag- 
in [0 a Gauss-BRST Ward identity was derived in the formal a = Coulomb gauge. 
This identity is the functional analog of the operator statement that the BRST symmetry 
transformation is generated by the Gauss-BRST charge. In the present section we shall 
show that this identity holds in the Coulomb gauge limit a — > from the Landau- Coulomb 
interpolating gauge. 

Consider the partition function in the Euclidean theory, with Coulomb type interpo- 
lating gauge, V • A + aA = (or = /). Having assured ourselves of Lorentz invariance, 
we set V = M = 0, and the partition function becomes 

Z(J, K,L) = J d$ exp[— £ + ($, J)], (78) 

where $ = (A, c, c, b), and 

£($, K, L) = S FP ($) + (K, sA) + (L, sc). (79) 

With E = / d 4 x A, the Lagrangian density reads 

A ee {l/A)Fl u + (K, + d'^c)D(A),c + L(-g/2) ■ (c x c) - %bA, + h\ (80) 

where d'^ = (ad)^ = (a<9 ,V). We are interested in the Coulomb gauge limit a — > 0. 
Because of the gauge arbitrariness of the Coulomb gauge discussed in the Introduction, 
this limit may be /3-dependent, with j3 = for the Landau- Coulomb gauge or /3 = a for 
the Feynman-Coulomb gauge. 

The Lagrangian density is BRST-closed, sA = 0. This implies the existence of an 
identity associated with the corresponding Noether current, which we now derive. For 
this purpose we make the infinitesimal change of variable of integration corresponding to 
a space-time dependent BRST transformation 

$^ = $ Q + e (x)s$ Q , (81) 
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where e(x) is space-time dependent, and a is an index the runs over all components of 
all integration variables. This change of variables leaves the measure rf$ invariant, and 
so, because e(x) is arbitrary, it yields the identity 

= J dQidtfn + sA^ + scJ c + bJ B ) exp[-S + ($, J)], (82) 

where is the Noether current of the BRST symmetry of A. If we integrate this identity 
over all space-time, the term d^j^ is annihilated, and we obtain the Zinn- Justin equation 
used previously. Instead we integrate over 3-space only, with spatially periodic boundary 
conditions, and obtain 

/ d " x{JA »W» + ^ " Je f| } = d oJ d< ^Q ex Ph S + J)\- (83) 
The conserved BRST charge Q is calculated from 

e = / " + (sc) w + ( s5 ' im ] - (84) 

where the fermionic derivatives are left derivatives, which gives 

Q = J d 3 x[-cD t F 0t - (K + ad c)(sc) + abD c\. (85) 

We wish to express the BRST charge in a way which will provide a Ward identity satisfied 
by the quantum effectve action T. For this purpose we observe that Q may be written 

Q = J d 3 x[-c^- + K ^-] + Q a (86) 

Q a = a J d 3 x s(bA - <9 cc) (87) 

is proportional to a, and is the integral of a BRST-exact density. 

The quantity jj- is the left-hand side of Gauss's law. In a canonical formulation, it 
is also the generator of local gauge transformations, so the first term of Q has the form 
of the generator of an infinitesimal gauge transformation with generator —c(x). For this 
reason, the last expression for the BRST charge Q remains correct if coupling to quarks 
is included in the Lagrangian density, and also in the phase-space representation which 
we shall introduce in the following section. 

From this expression for Q we obtain 

5Z 5Z 8Z S n r .« , , SZ „ 8Z S 



where 
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The expectation-value (Q a ) is calculated in the presence of all sources. In terms of the 
generator of connected correlation functions, W(J,K,L) = \nZ(J,K,L), this identity 
reads 

SW SW , 5W S „ r , , 5W „ 5W„ 



I 



dW OW OW f , d W OW 

d x{jA »W~ + Jc Jl ~ Jz 17 ] = d °J d x{jA ° JT " Ko ir ] + {Qa) - (89) 



We make the Legendre transformation to the quantum effective action K, L), which 
satisfies 

r j3 ,6T 5T 5T5T 7 5T. . r j3 . 5T ^ ST. ,■ , 

7 d + 5^ + = d °J dx{c jT - ^ - (90) 

Because Q a is proportional to a, one has Q a = in the formal Coulomb gauge a = 0. 
However Feynman integrals diverge in the limit a — ^ 0, so a precise evaluation is required 
to determine whether or not (Q a ) really vanishes in the limit a — > 0. In the following 
section we study this limit from the Landau-Coulomb interpolating gauge, (3 = 0, by 
means of a phase-space representation. By power counting of the k integrations, it 
is found that the correlation functions with dimensional regularization are finite in the 
limit a — > 0. It is found that, although (Q a ) does not in fact vanish linearly with a, 
nevertheless it does vanish like 

(Q a )=0{a 1 ' 2 ). (91) 

in the limit a — > 0. (See remark 3 at the end of the following section.) 

We now take the limit a — > 0, and set (Q a ) = 0. Only first functional derivatives of 
T appear, so the unacceptably singular expression of correlation functions at coincident 
points is absent, and this identity imposes a constraint on the renormalization constants 
of the elementary fields. As before, the Lagrangian multiplier fields b and c may be 
eliminated by means of their equations of motion, and the Gauss-BRST identity simplifies 
to 

, 3 ,6T 5f 5f5f n r j3 . ST T 5f^ 
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According to our results on renormalization, the quantum effective action f is finite 
when expressed in terms of renormalized quantities, 

f (X) = f r (X r ), (93) 

where X = (A,c,K,L,g,A) and X r = (A r ,c r , K r , L r , g r , fi). Here A is the usual ultra- 
violet regularization parameter, and \x is a renormalization mass. The renormalization 
constants satisfy 

Zk, = zi* z L = z;\ (94) 
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Moreover, for the Coulomb gauge, by rotational invariance, the matrix Z A is given by 
Z A ^ = di ag (Z Ao ,Z A ,Z x ,Z A ), and Z K ^ = diag(Z^\ Zz\ Zz\ Zz 1 ). Consequently the 
Gauss-BRST identity reads 

fd 3 x( ^ r ^ r I $r r 8r r ^_ Zc dy f d 3 x(c ^ r K o— ) (95) 
J 5A rfl 8K riX 5c r 5L r Z Ao J 5A r}0 ' 5L r 

Since all other quantities in this equation are finite, the ratio Z c /Za must also be finite. 
This implies that in the recursive renormalization procedure described above, the diver- 
gent parts of Z c and Za are equal in each order n. The iterative renormalization may 
be done so the finite parts are also equal in each order, and the equality 

Z Ao = Z c (96) 

is maintained. 

For this purpose we must show that the renormalized action E r also satisfies the 
Gauss-BRST identity. It is instructive to first verify directly that E satisfies this identity. 
Indeed by Noether's theorem the variation of E under the above space-time dependent 
BRST transformation is given by 

d i xe{x)d lx j lx (97) 
where is the Noether current. On the other hand we have 

<5E = [ d 4 xe(x)s^ l ^- 
J 8$i 

f l4 , 5E 5E 5E5E 7 5E X 
= / d 4 xe(x)( — — + — — + 6—). 98 
J 5K„ 5A„ SL 5c 5c 



Since e(x) is arbitrary, it follows that E satisfies, 

5E 5E 5E5E ,<5E 



SK^ + jLTc +b Tc = - d ^ <"> 

Upon integrating this equation over 3-space and using the above expression for the BRST 
charge Q, we obtain 



/ 



rt ,8E 5E 5E5E 7 5E N r . 5E tfE, . . , 

^^TTTir + + = / ^(c-- - ^ 077 ) - Qa- (100) 



•SApSKn 5c 5L 5c' .) y 5A n U 5L' 

4, 



We now introduce E(A, c, K,L) — J d x A, where 

A(A, c, X, L) = (1/4)F* + K^ l D(A) IJj c + L(-<?/2) • (c x c), (101) 
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so 

A(A, c, c, b, K, L) = A(A, c,K + d'c, L) - d'^bA^ + ^b 2 . (102) 
By the above reasoning we conclude that E satisfies the functional identity 



r , 3 ,52 5E <5E5E N . r , 3 , 5E 5E X , ino . 



5E 5E . <5E5E N 5E „ <5E. 

If one makes the change of variables 



A^ — Za^A t ^ 
L = Z c * L r 

9 = Z y g r 

E(A, K, c, L, g) = t r (A r , K r , c r , L r , g r ), (104) 

with Za q = Z c , this identity remains unchanged, so E r satisfies the same functional 
identity as E. This is the required condition for recursive renormalization. 

We have taken the limit a — > from the Landau-Coulomb interpolating gauge, (3 = 0, 
for which the estimates of the following section hold. In this gauge, as noted at the end 
of sect. (6), the renormalization constant c 3 = 0, so Z g Z c = 1. We therefore obtain in 
the a — > limit from the Landau-Coulomb interpolating gauge 

Z g Z Ao = 1. (105) 

Consequently the field gA is invariant under renormalization 

gA = svA-,o, (106) 

as are its correlation functions, including in particular the zero-zero component of the 
gluon propagator, 

D 00 (\x\,t) = g 2 {A (\x\,t)A (0,0)). (107) 

This quantity is independent of the cut-off A and the renormalization mass //, and conse- 
quently it can depend only on physical masses such as Aqcd- This holds for the instanta- 
neous part of D Q0 (\x\,t). However the instantaneous part of D 00 (\x\,t) may not be easy to 
separate uniquely (for example even in finite orders of perturbation theory), and a more 
accessible quantity is {7(|a?|) = — / dtD^(\x\, t). It also depends on physical masses only, 
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~ — * ~ — * ~ — * 

as does its fourier transform C/(|A;]) which is given simply by i7(|A;|) = Doo(k, &o)|fc =o- 
We write U(\k\) = g 2 c /k 2 . Here g c = g c (\k\/ A-qcd) is a running coupling constant defined 
in the Landau- Coulomb gauge that depends only on Aqcd- Such a quantity cannot be 
extracted from the gluon propagator in covariant gauges. Indeed to extract it in covariant 
gauges one must consider the Wilson loop which involves n-point functions of all order n. 



9 Coulomb gauge limit 

We now turn to a more precise analysis of the behaviour of the correlation functions 
when the Coulomb-gauge limit a — > is taken from the Landau- Coulomb interpolating 
gauge, characterized by (3 = 0. Because the gauge parameter a provides a rescaling of 
the time, instantaneous interactions appear as a approaches 0. 

Consider the partition function in the Euclidean theory, with Landau- Coulomb type 
interpolating gauge, V • A + aA = 0, 

Z = J d 4 Adcdcdb exp{ - J d 4 x [{^){E 2 + B 2 ) + ib(V ■ A + aA ) 

+ (atD c + Vc- Dc)}}, (108) 

where t = x represents Euclidean "time", = A\ — D^Aq] B = B(A), = D^(A). 
(The i appears in front of b, because b is here integrated over a real instead of imaginary 
contour.) For simplicity, we have suppressed all sources, and a summation on color 
indices is understood. 

We use the Gaussian identity exp[(-±) / d 4 x E 2 ] = J rf 3 Pexp[- / d 4 x (iP ■ E+(\)P% 
to obain the phase-space representation 

Z = J d*Ad 3 Pdcdcdb exp(-S), (109) 

— * — * 

where A and P are canonical variables, and 

S = J d 4 x[iP ■ {A- DA ) + (i)P 2 + (i)B 2 

+ib( V • A + aA ) + {aW c + Wc-Dc)}. (110) 

— * 

The phase-space action is BRST-invariant, with P transforming according to sP" = 

jabdpb^d 

We now make a linear change of field variable in order to diagonalize the gluon 
propagator, while keeping the action local. We pose 



A = VV (HI) 
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— * 

for which dA = const dip, and we shift A by 

A = A'-aVip, (112) 
for which d 3 A = d 3 A'. This simplifies the Lagrange-multiplier term 

ib(V ■ A + ad A ) = ib (V ■ A') (113) 

— * — » 

so it imposes the time-independent constraint V ■ A = 0, and we have 

iP ■ (A- DA ) = iP ■ (A> - aV^ - DA ), (114) 

— * — * — * — * — * —* • — * 

where D = D{A) = D{A' — aVip). We similarly separate P into its transverse and longi- 
tudinal parts, while keeping the action local, by introducing another lagrange multiplier 
field by means of the identity, 

const = JdQ5(V-P + V 2 n) 

= J dVldvexp[-i J d A xv (V • P + V 2 fi)], (115) 

which we insert into the partition function. We shift P' according to 

P = P'-Vtt, (116) 

under which d 3 P = d 3 P', so the new Lagrange-multiplier term becomes 

iv(V ■ P + V 2 Vl) =i v (V-P'), (117) 

and enforces the time-independent constraint V • P' = 0. The field Q represents the 
color-Coulomb potential. 

The partition function now reads 

Z = J d 3 A'd 3 P'dbdvdipdndcdc exp(-S") (118) 

where 

S' = [ d*x [iP ■ (A - DA ) + (i)P 2 + ()-)B 2 + ivV ■ P' + ibV ■ A' 

+(atD c + Wc- Dc)}, (119) 

B = B(A) = B(A' - aVip), and P = P' - Vfi. The first term in S' is given by 

iP ■ (A- DAq) = iP' ■ (A 1 - aV^ - DA ) 

-iVn-(A'-aV4>-DA ). (120) 
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To cancel cross terms in S' we shift the Lagrange multiplier fields, 

b = b' + tt 

v = v ' - ( a ^ + A ) -itt, (121) 

with dbdv = db'dv', and obtain, after integrating by parts in space and time and writing 
V 2 ^ = Ao, 

S' = J d A x [ iP ■ {A! -gAx A ) + i(aClA + Vfi • DA ) 
+ (\)p' 2 + (^)(Vfi) 2 + {\)B 2 

+iv' V • P' + ib'V ■ A' + (atD c + Vc-Dc)}. (122) 
The remainder of this section is an analysis of the action S' . The Lagrange multiplier 

— * — * — * — * 

fields b' and v' enforce the time-independent constraints V • A' = and V • P' = 0, on 

— * — * 

the canonically conjugate variables A' and P', and we call these "the transverse fields". 
The bose fields A and Q form a pair similar to the pair of fermi fields c and c, and we 
call this quartet "the scalar fields" . 
The corresponding free action 

S = J d A x [ iP' ■ A> + (i)P' 2 + (^(e^V^) 2 + iv' V • P' + i&' V • A' 
+ i (afii + Vfi • VA ) + (i)(VQ) 2 

+ {acc + Vc-Vc)]. (123) 

determines the free propagators. In momentum space the propagators of the transverse 
fields are given by 

D A >.A>. = (6 ij -k i k j )( k o + &)~ 1 
Dp, P . = (d^k 2 - k-kMkl + k 2 )- 1 

D P , A , = iko^ - kikj)^ 2 + k 2 )-\ (124) 

whereas the propagators of the scalar fields are given by 

D A() n = {ak 2 + k 2 r 1 
Dnn = 
D AoAo = k 2 {ak 2 + k 2 )- 2 

D c - c ={akl + k 2 )-\ (125) 
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Propagators of the ip field are obtained from ip = (V 2 ) _1 A . The new fields have conve- 
niently diagonalized the gluon propagator by separating the 3-dimensionally transverse 
and scalar parts. The transverse propagators have denominators (k 2 . + k 2 ), whereas the 
scalar propagators have denominators (a/c 2 , + k 2 ). Thus the scalar fields have a reaction 
time of order a 1 / 2 which is very rapid as a aproaches 0. Consequently it is natural to 
integrate out, if possible, the scalar fields and obtain an effective theory for the transverse 
degrees of freedom. 

To study the limit a — > 0, we separate the action into 4 terms, 

S' = S r + S X + S Y + S z . (126) 

Here S r is the free action So plus all vertices that are independent of a, 

S r = J d 4 x [ iP ■ (A> - gA' x A ) + (i)P' 2 + {\)B> 2 

+iv'V-P'+ ib'V-A' 

+i (aflA + ■ D'A ) + d)(Vfi) 2 

+ (afe + Vc- D'c)}, (127) 

— * — * — * — * — * — * 

where B' = B(A') and D' = D(A'). We shall see that S r has graphs that diverge as 
a — > 0, but that they cancel by virtue of an r-invariance. The term Sx consists of all 
vertices with 3 scalar fields and one power of a, 

S x = ag J rf 4 x[-iVjfi(V^ x A ) + t(A x c) - Vic(V^ x c)]. (128) 

The term Sy consists of a vertex with one power of a, 

S Y = agi J d 4 xP-(Viip x A )]. (129) 

There remains 

S z = f d*x(l/2)[B 2 (A' - aV^) - B 2 (A% (130) 

We first discuss the theory defined by S r , temporarily ignoring the vertices Sx, Sy 
and Sz that vanish with a. The action S r is at most quadratic in the scalar fields. Its 
vertices contain no powers of a and no time derivatives, so in momentum space there 
are no factors of ko at the vertices of S r . Consider a closed loop that consists entirely 
of scalar propagators with denominators (ak 2 + k 2 ). It is controlled by a time scale of 
order a 1 / 2 . The loop integral on k is effected by the change of variable ko = a~ 1//2 /c , 
which effectively eliminates a from the denominators, but the volume element of the loop 
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integral changes by dk = a~ l l 2 dk' . We conclude that each closed loop that consists 
entirely of scalar propagators and vertices of S r diverges like cr 1 / 2 . 

Nevertheless the theory defined by S r is finite as a — > 0, as we now show. We write 

S r = SV,i + SV,2 (131) 

where SV,i consists of all terms that contain only transverse fields and their Lagrange 
multipliers, 

Sr,l = Jd*x[ l P>.A'+ (i)P' 2 + (i)5' 2 

+iv'V-P'+ ib'V-A']. (132) 

It is independent of a. The remainder SV,2 also depends on the scalar fields and on a. 
It is helpful to express SV,2 in terms of the color charge density of the transverse fields 
p = gP[ x A' { , and the Faddeev-Popov operator M = — a<9g — V • D' characteristic of S r . 
We have, in an obvious notation, 

S r>2 = -i(p, A ) + MAq) + (1/2)(VQ, Vfi) + (c, Mc). (133) 

If one integrates out the ghost fields c and c, one obtains the Faddeev-Popov determinant 
det M. If one next integrates out A , one obtains S(MQ — p), which expresses the form 
of Gauss's law appropriate to S r . Finally the integral on dVt absorbs the Faddeev-Popov 
determinent 

J dVL&etM 5(Mn-p)exp[-(l/2)(Vn,Vf2)] = const, x exp(-S' coul ), (134) 

where 

S conl = (1/2)(VM"V, VM"V), (135) 

depends on the transverse fields only. It represents the non-local color-Coulomb inter- 
action, regularized however by the finite value of a. Thus the theory described by the 
local action S r = SV,i + 5V,2 that contains the scalar fields is equivalent to the theory with 
transverse degrees of freedom only, described by the non-local action SV,i + Sc OU i- 

Moreover SV,i + S cou \ at finite a provides a regularized version of the canonical action, 

The canonical action S can is a function of the canonical variables which are the transverse 
fields A' and P'. It is obtained by formal canonical quantization in the Coulomb gauge, in 
which one solves the constraints to eliminate the so-called unphysical degrees of freedom. 
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(136) 



To show that the theory described by the local action S r , or equivalently by SV,i+ Seoul, 
is finite in the limit a — > 0, we observe that the perturbative expansion of S cou i produces 
ladder graphs, in which the instantaneous parts are the horizontal rungs, corresponding 
to the instantaneous color-Coulomb interaction. Since these ladder graphs do not contain 
any instantaneous closed loops, they are finite in the limit a — > 0. To summarize: in the 
theory described by S r , each closed loop of bose and fermi scalars diverges like a -1 / 2 , but 
they precisely cancel to give a result that is finite as a — > 0. 

It is helpful to exhibit the cancellation between bosons and fermions in the theory 
described by S r by means of an r-symmetry. We express the action SV,2 in terms of the 
field Cl = Vl- M~*p, 

S r , 2 = MA Q ) + (c, Mc) + (1/2) (VO, VO) + (Vfi, VM"V) + Seoul- (137) 

Let r be a BRST-type transformation that acts on the scalar fields according to 

r Aq = c rc = 

rc=-itt rtt = 0, (138) 

— * — * 

and that annihilates the transverse fields and their Lagrange multipliers, rA' = rP' = 
rb' = rv' = 0. It is nil-potent, r 2 = 0. The action S r ^ may be written 

S r ,2 = Seoul + rtt, (139) 

where 

* = - (c, MAq) + (i/2)(Vc, Vfl') + i(Vc, VM _1 p), (140) 

and we have 

S r = 5 r)1 + + r*. (141) 

The first 2 terms depend on the transverse fields only, and are thus r-invariant, rSV,i = 
'"Seoul = 0. The last term r\l>, which contains all the dependence on the scalar fields, is 
r-exact. We have rS r = 0, and r is indeed a symmetry of the theory defined by S r . Now 
consider the integral over the scalar fields while the transverse fields and their Lagrange 
multipliers are held fixed. The effective action for the scalar fields is r^, which is r-exact. 
A theory whose action is exact under a BRST-type transformation is called "topological" , 
and has the property that the partition function, 

J dA Q dVLdcdc exp(-rW), (142) 
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is constant under continuous variations of the external parameters, namely the transverse 
fields and the parameter a. We have obtained the previous result, with the understanding 
that the cancellation of bose and fermi loops that diverge in the limit a — > is preserved 
by the r-symmetry of S r . The r-symmetry which transforms A into c explains the 
equality of renormalization constants that holds in the limit a — > 0, Za = Z c , which was 
established in the last section . 

We now come to the remaining vertices of S', namely Sx, Sy and Sz- These vertices 
formally vanish in the limit a — > 0, and they would not appear in formal canonical 
quantization in the Coulomb gauge. However because, as we have seen, there are closed 
loops in the expansion of Sz that are of order a -1 / 2 (and that cancel pairwise), we must 
verify whether insertions into these loops of the vertices Sx, Sy or Sz may give a finite 
result. These vertices are not r-invariant, so if there are such contributions there is no 
reason to expect that they cancel. 

Consider first the vertices of Sx which we call X-vertices. (Similarly we call r- 
vertices the vertices of S r etc.) The X-vertices are linear in a. They also contain one 
time derivative, so in momentum space they contain one power of ko = a _1 / 2 A; . Thus 
overall when an X-vertex is inserted into a closed loop of scalar propagators it gives a 
contribution of order a 1 / 2 . As we have observed from dko = a~ l l 2 dk' , the volume element 
for a closed loop consisting of scalar propagators is of order a -1 / 2 . Thus the presence of 
a single X-vertex in a closed loop of scalar propagators and r- vertices would give a finite 
limit, except for the fact that such a loop is odd in k' at large k' Q , and consequently a 
closed loop of scalar propagators with a single X-vertex is reduced to order a 1 / 2 . By the 
same reasoning, a closed loop of scalar propagators that contains two X-vertices (and is 
thus even in k' ) is also of order a 1//2 . Thus a single closed loop with one or two X-vertices 
vanishes like a 1 / 2 as a — > 0. However a closed scalar loop with two X-vertices has two 
external scalar lines, because each X-vertex is trilinear in the scalar fields. Consequently 
such a loop may be inserted into a closed scalar loop whose remaining vertices are all 
r- vertices. [See fig. (1).] This gives a two-loop graph, with two X-vertices, each of order 
a 1//2 , and two closed loops of scalar propagators, each of order gT 1 / 2 . This is finite in the 
limit a — > 0. (Further insertion of X-vertices gives a vanishing contribution in the limit.) 
We conclude that scalar bose or fermi closed loops do not decouple as a — > 0, but give 
a finite two-loop graph. This contribution is missing in formal canonical quantization in 
the Coulomb gauge. 

The analysis of the vertices of Sy is similar. Each F-vertex contains two scalar fields 

— * 

and one P' field. It also contains one power of a and one time derivative, so a F-vertex is 
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also of order a 1 / 2 . Again, insertion of single Y- vertex into a scalar closed loop would be 
finite except that it is odd in k' . We cannot connect up two Y-vertices by an additional 
scalar propagator because Y-vertices are bilinear in the scalar fields. However the Y- 
vertex contains the P' field which has the P[ — A'j propagator Pij(k)k (kl + k 2 )' 1 that 
contains k in the numerator. (It is the only propagator with k in the numerator.) Now 
consider a closed loop that consists of scalar propagators and one P' — A' propagator. All 
the vertices are r-vertices except for one Y-vertex at one end of the P[ — A'j propagator. 
[See fig. (2).] When the loop momentum ko is of order a~ 1//2 , the Y-vertex is of order 
a 1//2 , the P[ — A'- propagator is of order a 1 / 2 , and the volume element of the loop integral 
is of order a" 1 / 2 , so overall this closed loop is of order a 1 / 2 . However it has two scalar 
external lines that emerge from the two ends of the P' — A 1 propagator. Consequently 
this closed loop, which is of order a 1//2 , may be inserted into in a scalar closed loop 
consisting of S r vertices which is of order a -1 / 2 . [See fig. (3).] This again gives a finite 
two-loop contribution that is missing in canonical quantization in the formal Coulomb 
gauge. (Further insertions of Y-vertices give a vanishing contribution in the limit.) 

Finally, the vertices of Sz give vanishing contribution in the limit a — > 0, because 
when they contain 2 or 3 scalar fields they also contain 2 or 3 powers of a respectively. 

We summarize the results of this section: (1) The diagrams for which the ko integra- 
tions would diverge in the Coulomb-gauge limit, a — > 0, have been have been shown to 
cancel at finite a. The remaining diagrams are finite in this limit by power counting of 
the ko integration. (2) There are two-loop graphs of the scalar particles A -Q and c — c 
that are finite in the limit a — > 0, and that are missing from canonical quantization in 
the formal Coulomb gauge. It remains a logical possibility that these graphs are mere 
gauge artifacts that do not contribute to a gauge-invariant expectation-value such as a 
Wilson loop. However there is at the moment no argument to show that this is true. 

Remarks 

I. The correlation functions that do not involve the field P are the same as in the 
configuration-space representation, so the finiteness of the unrenormalized correlation 
functions in the Coulomb-gauge limit of the Landau-Coulomb interpolating gauge also 
holds in each order n for the configuration-space correlation functions. This implies 
that the configuration-space generating functionals Z, W , T and T are also finite in the 
limit a — > 0. Here an ultraviolet dimensional regulator e is understood to be in place. 
For the diagrams we have examined, the a-dependence at small a is given by a~ m ^ 2 , 
where m is a no n- negative integer. (The terms with negative powers of a 1 / 2 cancel.) 
These powers are e-independent, and so cause no trouble in the e — > limit (as would, 
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for example, terms like a e ). This is because the terms that diverge with a come from 
divergences in the one-dimensional ko integrations and are not affected by dimensional 
regularization which is a continuation in the number of spatial dimensions. Likewise 
the cancellation of terms that diverge as a — > is assured by r-invariance, and is also 
dimension-independent. Moreover rg iv (a, e), eq. (|67|) , has a simple pole structure in e. 
Consequently the finiteness of T n (a,e) = r^(a, e) + r^ iv (a,e) as a — > implies that the 
residue of T^ iv (a, e ) an d r^(a, e) are separately finite as a — ► 0. Although we have not 
made an exhaustive examination of all diagrams, we expect that the remaining diagrams 
behave similarly, and thus that the renormalized correlation functions are finite in the 
Coulomb-gauge limit of the interpolating Landau-Coulomb gauge. 

2. We may regard the finite value of the 2- loop scalar graphs that are missing in the 
formal Coulomb gauge, a = 0, as an anomaly of the r-symmetry; for the action S'(a) is r- 
invariant at a = 0, rS'(0) = 0, but not at finite a, and the symmetry is not regained in the 
limit a — > 0. This comes about because individual graphs diverge in this limit, and they 
combine with subgraphs containing r-noninvariant vertices which vanish in the limit, to 
give a finite result. However the divergent graphs result from a part of the action r^(a) 
that is r-exact at finite a, and thus topological. This assures that the divergent graphs 
cancel each other, so that the limit is finite. It also preserves the equality, Za = Z c , 
among the limiting renormalization constants found in the last section, eq. fl96|), which 
would hold if the transformation r A = c were actually a symmetry of the limiting theory. 

3. To establish the Gauss-BRST Ward identity of the last section, there remains to 
verify that (Q a ) = in the limit a — > 0, where the expectation value is calculated in 
the presence of all sources. Here Q a , is the part of the total BRST charge Q defined in 
eq. fl87j), 

Q a = aJ d 3 x[bD c - d bc + d c{-g/2)c x c]. (143) 

To evaluate (Q a ), one makes a diagrammatic expansion of each term by the method of 
the present section, using b = b' + Q. The only non-zero propagator of the b' field is the 
b' — A[ propagator, k^ak^ + k 2 )~ l . For example, consider the contribution of the term 
b c d = b' d c d + Cl d c d . The term Q d c d looks dangerous because it contains two powers of 
ko- However the only non-zero propagator of the Q field is the Q — Aq propagator, and 
the vertex where A$ is absorbed is proportional to a. A typical graph representing the 
contribution of Q d c d to the fourier transform / dtexp(ipot)(Q a (t)) is illustrated in fig 4. 
(Note that Q a is a color scalar, so it must decay into at least two quanta, namely c and 
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Ai in fig. 4.) This graph contributes 

ap (2vr)- 4 y d A k k 2 (ak 2 + k 2 )~ l a(k + q )[a(k + q ) 2 + (k + q) 2 }- 1 

xkMh+Po^ + k 2 ]- 1 . (144) 

We rescale the variable of integration dk = a^^dk^ and obtain a contribution of leading 
order a 1 / 2 , keeping in mind that terms that are asymptotically odd in k' are suppressed 
by a 1 / 2 . The other terms are evaluated similarly. One finds that each term of Q a gives a 
contribution to the expectation- value of order a 1 / 2 . QED. 

4. We we have seen that a closed loop of unphysical particles, with propagators 
(k\ + k\ + k\ — a/cp) -1 , is of order a -1 / 2 as the Coulomb-gauge limit of the Landau- 
Coulomb interpolating gauge is approached. This behavior came from the rescaling, 
k = a~ l / 2 k' 0) that makes the loop integral of order dk = a~ l ^ 2 dk' . However in the light- 
front or axial gauge, the unphysical propagator is [a(k\ + k 2 ,) + (1 + a)(k 2 — /cq)/2] _1 or 
[a(/c 2 + /c| — kfy + k 2 }' 1 , and the required rescaling, (k 1: k 2 ) = a~ 1 ^ 2 (k' 1: k' 2 ) or (k , ki, k 2 ) = 
a~ 1 / 2 {k' Q , k[, k' 2 ), gives an uphysical closed loop integral of order a -1 or a~ 3 / 2 . Thus the 
light-cone and axial gauge limits appear to be more singular than the Coulomb-gauge 
limit, and additional cancellations would be required to give finite correlation functions. 

10 Conclusion 

We briefly review our results. We have addressed the problem of the existence of "phys- 
ical gauges", by the device of interpolating gauges which interpolate linearly between 
a covariant gauge, such as the Feynman or Landau gauge and a physical gauge such 
as the Coulomb or light-cone gauge. For example, the interpolating Landau- Coulomb 
interpolating gauge is defined by the gauge condition ad A + V • A = 0, which gives the 
Landau gauge for a = 1, and the Coulomb gauge is achieved in the singular limit a — > 0. 
More generally an interpolating gauge is defined by the condition a^ u d u A^ — (or = /), 
where a is a non-singular numerical matrix, and a "physical" gauge is a limiting case in 
which a becomes singular. 

In general the interpolating gauge breaks Lorentz invariance as well as local gauge 
invariance. Nevertheless we are able to establish the existence of the perturbative expan- 
sion and perturbative renormalizability of the interpolating gauges in full generality, by 
extending the BRST method to include the Lorentz group in addition to the usual local 
gauge group. This extension is necessary to control the form of divergences, for example 
to show that the divergent coefficients of the term ceE 2 + cb§ 2 are equal, ce = cb- 
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The enumeration of the possible divergence terms that are BRST-invariant is not sub- 
stantially more difficult than for Lorentz-covariant gauges. Moreover the matrix a is 
a gauge-parameter in the sense that the expectation values of physical observables are 
independent of a, as long as a is non-singular. Thus the interpolating gauges are strictly 
equivalent to the covariant gauges. 

However the singular limit to a physical gauge is quite subtle. It is analyzed in the 
present article for the Coulomb gauge limit, a — * 0, from the Landau- Coulomb interpo- 
lating gauge. There are closed bose and fermi-ghost loops that become instantaneous 
in the limit a — > and that individually diverge like a -1 / 2 . We we use a phase space 
representation and a linear shift of field variables to exhibit the cancellation of loops that 
diverge like a _1//2 , and to show by power counting of the fco integrals that this limit gives 
finite correlation functions. 

An important aspect of this limit is that there are also closed bose and fermi one-loop 
graphs that are not present in the formal Coulomb gauge (a = 0). Although they vanish 
like a 1 / 2 , they cannot be neglected because, when these one-loop graphs are inserted into 
the above-mentioned closed loops that diverge like a -1 / 2 , they give a finite contribution. 
Consequently the closed bose and fermi-ghost loops do not decouple in the Landau- 
Coulomb gauge limit, but give a finite two-loop contribution. 

One logical possibility is that these two-loop ghost contributions are merely a gauge 
artifact that do not actually contribute to expectation-values of gauge-invariant objects 
such as Wilson loops. However there is at present no argument in hand to show this. If 
these two-loop bose and fermi-ghost graphs do contribute to physical expectation values, 
then the traditional picture of the Coulomb gauge would have to be revised. The state 
space would not be simply describable in terms of transverse gluons. In the latter case 
the Coulomb gauge is not more unitary than other gauges, in the sense that it cannot 
be simply described in terms of the classical dynamical variables that remain after the 
constraints are solved. Indeed we are unable to provide a set of Feynman rules to be used 
in the Coulomb gauge at a = 0, although we have shown that both the unrenormalized 
and renormalized correlation functions are finite in the limit a — > of the Landau- 
Coulomb interpolating gauge. 

Nevertheless there is a reward to be gained by taking this limit. For we have shown 
that the Gauss-BRST Ward identity holds in the Coulomb gauge limit of the Landau- 
Coulomb interpolating gauge. This identity is the functional analog of the operator 
statement that the BRST symmetry transformation is generated by the Gauss-BRST 
charge. Among other things, it implies that gA Q is invariant under renormalization, 
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gAo = g^Aft. This means that all correlation functions of gAo are renormalization- 
group invariants, including in particular the time-time component of the gluon propagator 
g 2 D 00 . It depends only on physical masses such as Aqcd, but is independent of the cut- 
off or the renormalizaion mass. Thus the Coulomb-gauge limit of the Landau- Coulomb 
gauge provides direct access to renormalizat ion-group invariant quantities, whereas no 
component of the gluon propagator has this property in covariant gauges. Indeed in co- 
variant gauges one must go to the Wilson loop, which involves gluon correlation functions 
of all orders, to obtain a renormalizat ion-group invariant quantity. For this reason the 
Coulomb gauge may prove advantageous for non-perturbative formulations. In partic- 
ular, the instantaneous part of g 2 Doo may be a confining color-Coulomb potential that 
may serve as an order parameter for confinement of color ||. 
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11 Figure Captions 

1. Diagram with 2 scalar loops and 2 X-vertices, and any number of r- vertices. 

2. One- loop diagram with a single P' — A 1 propagator and an r- vertex. 

3. Insertion of graph of fig. 2 into a closed scalar loop of r- vertices. 

4. A typical graph contributing to (Q a ). 
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